Angular Dependence of Jet Quenching Indicates 
Its Strong Enhancement Near the QCD Phase Transition 
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We study dependence of jet quenching on matter density, using "tomography" of the fireball 
provided by RHIC data on azimuthal anisotropy V2 of high p t hadron yield at different centralities. 
Slicing the fireball into shells with constant (entropy) density, we derive a "layer-wise geometrical 
limit" v™ ax which is indeed above the data V2 < v 2 nax . Interestingly, the limit is reached only if 
quenching is dominated by shells with the entropy density exactly in the near-T c region. We show 
two models that simultaneously describe the high p t «2 and Raa data and conclude that such a 
description can be achieved only if the jet quenching is few times stronger in the near-T c region 
relative to QGP at T > T c . One possible reason for that may be recent indications that the near-T c 
region is a magnetic plasma of relatively light color- magnetic monopoles. 

PACS numbers: 25.75.-q , 12.38.Mh 



Introduction. — Recent experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) are dedicated to study 
possible new forms of QCD matter, with increasing 
energy density. In such collisions the produced mat- 
ter equilibrates as Quark-Gluon Plasma (QGP)[l[ and 
then cools down through the near-T c (M) phase (M 
for mixed/median/magnetic Q) into the usual hadronic 
phase (H). To probe the created matter in an externally 
controllable way, like using X-ray for medical diagno- 
sis is impossible. However, high energy jets are inter- 
nal probes: propagating through the fireball, they in- 
teract - and thus obtain important information about 
the medium - as proposed long ago in Refs 0,1,(1. In 
heavy ion collisions this energy loss can be manifested 
in the suppression of observed hadron spectra at high 
transverse momenta pt, as well as in the suppression of 
back-to-back di-hadron correlations with a high-^ trig- 
ger, when compared with pp and d-A collisions. The "jet 
quenching" phenomenon is one of the major discoveries 
by the RHIC experimental program^. 

The suppression is quantified by comparison of the 
inclusive spectra d 2 N AA /dptd-q in ion-ion(AA) collision 
to a nucleon-nucleon(pp) reference d 2 a NN /dp t dr/ via the 
Nuclear Modification Factor Raa{pi)- ■ 



RAA{pt) = 



d 2 N AA /dp t dn 



T AA ■ d 2 a NN /d Pt d V 



(1) 



with Taa the nuclear overlap function which scales up 
single NN cross section to AA according to expected 
number of binary NN collisions without modification. 
Thus a Raa smaller(larger) than unity means suppres- 
sion (enhancement) due to medium effect. At RHIC this 
ratio at large p t > QGeV has been measured to be a 
constant, about 0.2 for the most central AuAu collisions. 
Accurate calibration of hard processes in pp and dAu 
collisions, as well as with hard photon measurements 
(which show no quenching) Q resulted in quite accu- 



rate knowledge of jet production geometry, for any im- 
pact parameter b (or centrality bins, often characterized 
by the number of nuclcon participants N part in a col- 
lision event). While quenching is firmly established as 
a final state effect, many efforts to understand its mi- 
croscopic mechanism are not yet conclusive. Those in- 
clude pQCD gluon radiation with Landau- Pomeranchuk- 
Migdal (LPM) effect 0], synchrotron- like radiation on 
coherent fields elastic scattering loss etc. The 
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fate of deposited energy was discussed in Refs 
led to predictions of "conical flow" correlated with ex- 
perimentally observed conical structures in correlations 
involving 2 or 3 particles, for reviews see e.g. (l3. 14|. 

Jet tomography and the geometric limit. — In non- 
central collisions the overlap region of two colliding nuclei 
has almond-like shape: thus jets penetrating the fireball 
in different directions lose different amount of energy ac- 
cording to their varying paths. Their yield distribution 
d 2 N / dp t d(f) in azimuthal angle (f> (with respect to the re- 
action plane) for high p t h adrons thus provides a "tomog- 
raphy" of the fireball 15, [l|| [l7j]. We will focus on the 
second Fourier coefficient 



d<P cos(20) [d 2 N/d Pt d(f>} 
Jo n d<f>[<PN/dp t dcl>] 



(2) 



depending on impact parameter b for large pt > 6 GeV 
where hard processes dominate and dependence on p t is 
weak 18]. 

Unexpectedly, measured t>2 {pt , b) happen to be consid- 
erably larger than what jet quenching models predicted. 
The aim of our work is to provide simultaneous descrip- 
tion of both Raa and vi at high p t based on theoreti- 
cally known geometry of jet production and bulk matter 
evolution. One important concept of the analysis is the 
so called geometric limit, first suggested by one of us in 
171 ]: the observed asymmetry should be less than some 
value V2(largept,b) < v™ ax (b) provided by the geome- 
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try of the overlap region of two colliding nuclei. The 
idea (l7j was that for very strong quenching only jets 
emitted from the surface of the almond can be observed. 
Two other assumptions were made, namely: (i) quench- 
ing is proportional to matter density; (ii) colliding nuclei 
were approximated by homogeneous sharp-edge spheres. 
However even early experimental data showed that V2 
is actually well above this bound. Subsequent studies by 
Drees, et al [19( relaxed the second assumption, with real- 
istic nuclear shapes, which only made contradiction with 
data even stronger (see e.g. their Fig. 3(d)). 

The main lesson from those studies is that quenching 
is not proportional to the matter density, but a nontrivial 
function of it. Assuming some form of this function, one 
can then calculate both observables vi (b) and Raa ■ 

Layer-wise geometrical limit. — Systematically slicing 
the (expanding) fireball into shells with the entropy den- 
sity s a < s < Sb, we calculate what B,AA{b) and 112(b) 
would result with such single shell being the sole source 
of quenching by a Glauber simulation of AuAu collisions 
and jet production as in 17, With the quenching 



function k(s) assumed to be concentrated at this slice 
K a b-9(s — s a )-6(sb — s), the distribution in survival proba- 
bility / can be calculated and directly leads to evaluation 
of Raa- 
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, Raa =< f n ~ 2 > 



a 8.10 

(3) 

Extra I in the path corresponds to radiative LPM the- 
ory • The power index n comes from the ttq v* spec- 
trum in pp collisions, see detailed discussions in [181]. For 
each density shell the absorption coefficient n a b (in unit 
fm) is then fixed by Raa data (ljj parameterized by 
Raa (pt > 5GeV) = [1 - 8.3 • 10" 3 • A^]"" 2 . Then 
we calculate V2, by sampling half of the jets travelling in 
x directions ±5° and the other half in y direction and 



AA 



m. For 



extracting the difference in the respective R 
the Glauber initial condition we follow hydro calculations 
(see e.g. (2(j) to scale entropy density with local partici- 
pant density, and for bulk evolution we use 1-D Bjorken 
dilution which is appropriate till time ~ 10/to/c (see 
e.g. [2l|). Jet production points are simulated according 
to binary collision density. We have 24 entropy shells, 
(0,1],(1,2],...,(23,24] (in /fm 3 units). 

The resulting V2 for three impact parameters b = 
5, 7, 10 fm (bottom-to-top) are shown in FigfJTa). 

(i) Note that certain entropy shells produce V2 much 
larger than the old geometric limits of refs.[I3> 03, cor- 
responding to surface emission (small s at the left side of 
the plot). 

(ii) the existence of the maximum u™ aa: (6) leads to layer- 
wise geometrical limit: its dependence on centrality is 
shown in FigQJb) by filled big blue diamonds. 

(iii) Interestingly enough, the entropy shells where the 
maxima occur (for all centralities) correspond to the 
same interval s = 4 — 8 /to -3 , which is in fact quite spe- 



cial: it corresponds exactly to the vicinity of the QCD 
phase transition (see e.g. [24|). These curves reflect not 
only the geometry of the respective entropy shells, but 
also their placement relative to the jet production points. 

After these studies of single shells, we turn to the com- 
piled high-p t RHIC data on V2(b), shown in FigJTJb). We 
include only data for "hard" hadrons with p t > 6 GeV 
from PHENIX (open green boxes) and STAR (open ma- 
genta stars) collaborations. Comparing these data points 
to our layer-wise geometric limit (filled big blue dia- 
monds), we do observe that all the data points are (within 
error bars) indeed below this proposed bound. We also 
show V2(b) lines which would come out if all jet quench- 
ing would be due to two other single entropy shells, with 
s = (ll,12]/m -3 (filled small purple diamonds) and 
s = (23,24]/m -3 (open blue diamonds). Those corre- 
spond to the QGP phase, near and far from the tran- 
sition region: the values of 1/2(6) from those shells are 
significantly smaller than the maximal. Now we qualita- 
tively understand the experimental trend: going from the 
more central to the more peripheral collisions, quench- 
ing geometry shifts from quenching at high density shells 
(QGP), to the near-T c region at N p ~ 100 (approaching 
the upper limit). For extremely peripheral collisions we 
expect V2 to decrease again, reflecting geometry of the 
low entropy density shells (the hadronic phase). 

Modelling tomography of jet quenching. — We now 
turn from individual shells to realistic models, describ- 
ing the combined effect of all of them. 
Model A — a two-phase scenario model, in which we 
assume the quenching function n(s) with two parame- 
ters: one in the near-T c region and the other for the 
QGP phase, i.e. 

K (s) = k r x [1 ■ 9(s - si) ■ 6{s c 2 -s)+X-0(s- s r 2 )] (4) 

with s\ = 3//m 3 and s\ — 11/ /to 3 bracketing the near- 
Tc region. The parameter kr is globally fitted from 
RAA{N P art) (for each given A), while A characterizes the 
relative quenching strength between the near-T c region 
and the QGP, with its best value to be determined from 
a global fitting for V2{N part ). 

Model B — a scenario featuring peaked quenching 
strength at T c , which assumes 



k(s) = k r x [e 



-(- 



a 



(8 -8l)+£ -0(8-81)] (5) 



with s c = 1 1 fm? and s c w = 2//to 3 spanning the near-T c 
region according to lattice results [241]. 

Schematic sketches of the two models' k are shown in 
FigHfleft) and \ 2 /D.o.F from fitting the v 2 data (both 
the PHENIX and the STAR points), with a variety of 
choices of A (Model A) / £ (Model B), are shown in 
FiglUnght). The plots suggest that current V2 data fa- 
vors the relative quenching strength A = 0.4 for Model 
A and £ = 0.2 for Model B, both favoring a scenario that 
in relativistic heavy ion collisions the jets are quenched 
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FIG. 1: (color online) (a) The V2 obtained for each entropy shell at b = 5/m(dashed),7/m(solid), and 10/m(dotted) res pec tively; 
(b) v 2 rlax for high pr hadrons calculated at different N par t as compared with available RHIC data from [181. |22| and [231] . 
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FIG. 2: (color online) (left) Schematic demonstration of the 
quenching functions of our Model A and B; (right) The 
X 2 /D.o.F when fitting the V2 data with different values of 
parameters A (£) in our Model A (B), see text. 



about 2-5 times stronger in the near-T c region than the 
higher-T QGP phase. 

We also plot in Fig[3] the V2(N part ) obtained with 
the above optimal parameters: Model A with kr — 
0.00435/m and A = 0.4, Model B with k r = 0.00745/m 
and £ = 0.2. Both of them describe current data very well 
and predict rapid dropping of v-i at the very peripheral 
end N p < 100. 

Conclusions and discussion. — We started with the 
calculation of the "layer- wise geometric limit" for models 
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FIG. 3: (color online) Comparison between V2 experimental 
data and V2 calculated from our models, see text. 



describing jet quenching 

v 2 (b)<vr x (b) 



(6) 



where the r.h.s. is shown by the filled big blue diamonds 
in Fig[ljb) and corresponds to particular density shells. 
Unlike previously proposed bounds, this one is indeed 
satisfied by all available data (for large enough pt , within 
error bars). The limit can be reached only when the jet 
quenching is overwhelmingly dominated by the matter 
shells with the entropy density s — 4 — 8 /m -3 since 
only those have the right geometrical properties: the data 
points suggest this seems indeed to be the case for the 
AuAu collisions at RHIC at ~ 100. 

While previous models [13, L3| failed to reproduce the 
high p t t>2 and Raa simultaneously, we now presented 
two models which can do so. The key is the nontrivial 
dependence of quenching on the (entropy) density. We 



4 



concluded that the angular dependence of jet quenching 
indicates its strong enhancement near the QCD phase 
transition, about several times stronger than in the QGP. 

Why can it be so? Perhaps a near-T c peak in jet 
quenching should not be too surprising, as we already saw 
similar peaks/sharp- valleys around T c for other proper- 
ties of QGP, from trace anomaly, specific heat and speed 
of soundjUl to shear and bulk viscosities [25]. Recently 
the jet quenching strength was found to inversely related 
to shear viscosity in weakly coupled QGP[26[ — such 
relation if naively extrapolated and combined with the 
minimum of shear viscosity at T c would also point to a 
near-T c peak of jet quenching. It was also proposed in 
27j that switching on quenching only after a global time 
T q rsj 2 fm one can obtain better values of the asymme- 
try: such effect is incorporated by near-Tc dominance in 
a much more plausible manner via local density evolution. 

A microscopic explanation may be provided by recent 
magnetic scenario for the near-T c QCD plasma, in which 
this narrow T-region is treated as a magnetic plasma of 
light monopoles^2j]. In the same region quarks/gluons are 
few times heavier and thus get less energy for the same 
momentum transfer. When a fast electric charge (the jet) 
penetrates such plasma, its strong transverse magnetic 
field easily accelerates the abundant light monopoles into 
an overheated magnetic "coil" behind it via the dual 
Faraday effect, leading to substantial energy loss of the 

jet m. 

It will be interesting to extend the present study to dif- 
ferent colliding nuclei A and beam energy y/s: the data 
are becoming available (see e.g. [2{|) and the phenomena 
are rich as the jet production, the bulk evolution, and 
the pp reference all scale differently with A and i/s. In 
Fig|3] we have included the prediction for high pt vi of 
CuCu 200GeV collisions from our model B fixed by AuAu 
(orange filled triangle) , to be tested by data. More dedi- 
cated studies (including different initial scaling, different 
path length dependence, etc) will be reported in [301 ]. 

Note added: After the Letter was submitted, 
PHENIX run7 preliminary data were released [22]. They 
are now included in FigfIJa) and Fig[3] (squares): as one 
can see they agree with our model well. As also shown 
in [22| , most other models of quenching give t>2 2-3 times 
smaller than data. 
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